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INTRODUCTION 


Tue sky component due to a given unobstructed window opening, at any 
point inside a room, in any plane, is defined as the ratio of the inteusity of 
daylight illumination, in that plane, due to the patch of sky visible from that 
point through the opening, to the simultaneous intensity in a horizontal 
plane at an external point open to the entire sky vault. This component is 
independent of the absolute brightness of the sky, but is dependent on its 
surface distribution. 


The general practice for design purposes is to take the worst conditions: 
and a heavily overcast sky defines the most adverse condition for daylighting. 
Prior to 1955, a heavily overcast sky was generally assumed to have a uni- 
form distribution of illumination over its surface, and was accordingly called 
a ‘uniform sky’.. But measurements carried out by Moon and Spencer (1942) 
indicated that for a sky entirely covered over with thick dark clouds, the inten- 
sity of illumination at any point of the sky was independent of the azimuth 
and the position of the sun at the time, and varied only with the altitude in 
accordance with an empirical approximate formula 


B, ogy Ut — 
where B, and B, are the intensities respectively at the zenith and at an alti- 
tude 6. Subsequently, the International Commission on Illumination (C.L.E.) 
in its sessional meeting at Zurich (1955) defined and recommended for general 
adoption a “‘C.I.E. Standard Overcast Sky” in which the intensity varied 
in strict accordance with the Moon-Spencer formula. Since then, the day- 
light illumination inside a room is expressed in terms of the intensities, both 
inside and outside, arising from this Standard Sky only. 


As early as 1924 Yamauti derived equations to give the sky components 
with a “Uniform sky” for vertical rectangular openings. A geometrical 
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method was employed to obtain the components at any interior point, in 
three planes, (a) horizontal, (b) vertical and perpendicular to the window 
plane and (c) vertical and parallel to the window plane. A need has arisen, 
since 1955, to derive similar equations for a C.I.E. Standard Overcast Sky. 
In the absence of such equations, Waldram’s diagram (1950), originally con- 
structed for a uniform sky as early as 1923, was modified by Walsh (1959) 
to suit the standard sky. Rivero (1958) used Yamauti’s equations in draw- 
ing up his tables for a uniform sky and later modified those values to give 
tables of sky components for a non-uniform standard sky. Hopkinson, 
Longmore and Murray Graham (1958) prepared certain Simplified Daylight 
Tables for glazed openings with a standard sky. But in both an indirect 
method of computation alone was used to allow for the non-uniform distri- 
bution of intensity. 


The present paper gives exact trigonometrical equations for direct cal- 
culation of sky components due to rectangular windows, both vertical and 
horizontal, for the C.I.E. Standard Overcast Sky. The same equations will 
enable direct calculation of sky components due to any rectangular sloping 
window also, as it can be resolved into a horizontal and a vertical opening. 
The method adopted is quite general. Yamauti’s equations can be derived 
by the same method. It could also be employed to derive sky component 
equations arising from an intensity distribution other than that of the 
standard sky, provided that intensity is a known function of the altitude and 
independent of the azimuth or vice-versa. 


GENERAL METHOD OF DERIVATION 


In Fig. | ABCD is the rectangular vertical window of height A and 
length /. P is a point on the normal to the window through A and distant d 
from it. A sphere is described with P as centre and radius d. A projection 
of ABCD on the spherical surface is AB’C’D’ which is proportional to the 
patch of the sky visible from P through the window. The intensity at P due 
that portion of the sky is the same as from AB’C’D’ if the same luminous 
intensities are assumed at corresponding points. 


The values of 8, 8’, y, y’ the angles subtended at P by the sides AB, DC, 
AD, BC respectively are given by the following equations: 


(a) tan 8 =! 


h 
(b) tan y = 7 
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(c) tan Bp’ = tan B cosy 
(d) tan y’ = tan y cos B (1) 


The following additional relationships between them can also be easily 
derived : , 


(a) sin B’ = sin 8 cos y’ 
(b) sin y’ = sin y cos fp’ 
(c) cos B’ cos y = cosy’ cos B (2) 


In Fig. 1, M is a point at azimuth y and altitude 6. N,, Nz, Nz are the 
feet of perpendiculars from M on the three mutually perpendicular planes, 
PAY, PZA, and PZY, respectively. A small element of area of the spherical 
surface at M between azimuths % and % + dy and altitudes @ and @ + d@ 
is considered. The area of the element is' d? cos 6 dé dj. If B is the lumi- 
nous intensity at M, the intensity of illumination at P due to the element, 
on a plane perpendicular to PM, is Bcos@d@d. This expression when 
multiplied by sin @ gives the component of the intensity in the horizontal 
plane PAY, when multiplied by cos N,MP = sin MPN, = sin ¢ cos @ (vide 
equation 2 a) gives the component on the plane PZA, and when multiplied 
by cos N;MP = sin MPN; = cos % cos 8 gives the component on the vertical 
plane PZY. B is a constant for the uniform sky and 


B= 3, = 3,1 +2804 





for the standard sky. The above expressions are integrated over the spheri- 
cal area AB’C’D’ first with respect to @ and then with respect to %, between 
appropriate limits, and later expressed as ratios of the corresponding hori- 
zontal intensity due to the total unobstructed hemisphere of the sky, to give 
the respective sky components. 


HORIZONTAL INTENSITIES DUE TO THE ENTIRE SKY VAULT 


Total horizontal intensity (Ty) at P due to the uniform sky hemisphere 
is given by 
Ti2 27 


Tu =Bfs f sin@cos@ dé dy 
0 


0 
= 7B (on integration) (3) 


Total horizontal intensity (Ts) at P due to the standard sky hemisphere is 
given by 
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viz 27 


= Be ff 1 + 2sin 0) sin 0 cos 0 0 ap 


Ts 3 


; 7Bz (on integration) (4) 


SKY COMPONENTS FOR VERTICAL RECTANGULAR WINDOWS 
Sky component (fy.q) at P on a horizontal plane with a uniform sky 


is given by 


0 B 
B f J sin @cos 0 d0 ds 
fia = 2° a vices 
u 





[Tu is given by (3)] 


where 4’ (ref. Fig. 1) is given by tan 6’ = tanycos% (vide Equation | d) 
where tan y is a constant. The above equation when integrated over the 
area AB’C’D’ and simplified gives 


e.38 B — B’ cosy (5) 


V.H In 














Sky component (Fy.,) at P on a horizontal plane with a standard sky 
is given by 








m— ero - 2 tS 











Equations of Sky Components with a “‘ C.I.E. Standard Overcast Sky” 237 
o sR 
= f f (1 + 2sin @) sin @ cos 6 dé dys 
0 


0 


Fys = T; [Ts is given by (4)] 





which on integration and simplification gives 


Fat = oc (B — B’ cos 7) + 2. sin“ (sin B sin y) 


‘a 7 (sin 2y sin f’). 6 


Similarly expressions for fy,, and Fy,, the sky components on a ver- 
tical plane through P perpendicular to the window, for uniform and standard 
skies respectively, and corresponding expression for fy,, and Fy,, the sky 
components on a vertical plane through P parallel to the window are derived 
and given below. 





fa =Y— cor (7) 


Fy. = 2 (y — ¥' cos B) + = (1 — cos B — cosy + cos B cos y’) 


(8) 
ae B’ siny as y’ sin B (9) 





se rae 2 #3 ow 
Fy.y = a. sin y + y’ sin 8) + Tn (sin 8B — sin B’ cos? y) (10) 
Equations 5, 7 and 9 have been given by Yamauti (1924). 


Sky COMPONENTS FOR HORIZONTAL RECTANGULAR WINDOWS 


Figure 2 gives the diagram of a horizontal rectangular opening ABCD. 
AB'C’D’ is its projection on the surface of a sphere described with P as 
centre and radius d, and is proportional to the patch of the sky visible through 
the window from P. 8, f’, y, y’ are again taken as the angles subtended 
by the sides AB, DC, AD and BC respectively. Their values will be given 
by the same equations (1) and (2). The value of the angle a (ref. to Fig. 2) 
is given by the equation 

tan y 
tana = fan B (11) 
For any element of area at M included between azimuths % and % + dy 
and altitudes 9 and + d@ the expression By cos @ dé di where By is a 


















constant for the uniform sky or equal to Bz (1 + 2 sin 6)/3 for the standard 
sky gives the intensity at P in a plane normal to the line joining P with the 
element M. This expression is multiplied by the factors sin 0, sin ¢ cos 8, 
or cos % cos 8 as before to give the components of intensity at P in the hori- 
zontal plane, 8 plane (APX), and y plane (APY) respectively. The inte- 
grations however are carried out in two stages, between limits 7/2 and @’ 
for 6, and a and 0 for % in the first stage, and between limits 7/2 and 6” for 
6, and 7/2 and a for % in the second stage, where 0’ and @” are functions of 
% given by the equations 


tan 0’ = cot Bcos¢ (12) 
tan 0” = coty sin ¢ (13) 


The two integrated values are added together and expressed as a ratio of 
the horizontal intensity due to the total hemisphere of the sky Ty for the 
uniform sky or Ts; for the standard sky (Equations 3 and 4) to give the res- 
pective sky components. The following equations can then be derived, 
where fun, fa.g and fy., are the sky components with an uniform sky and 
Fun» Fug and Fy., are the sky components with a standard overcast sky, in 
the horizontal, 8 and y planes respectively. 


With uniform sky 


fast = 7 (B' sin y + 7’ sin B) (14) 
fae = >, (B — B' cosy) (19 


facy = 30 (9 — ¥ 008 B). 
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With C.ILE. Standard overcast sky 


Fun = 7ge (6 sin y + y’ sin 8) + 5 + 2 (sin 28 sin y 


+ sin 2y sin p’) — rd sin“ (sin a cos f) 


2 
— sin! (cos a cos y) (17) 
Fue = yp (B — B’ cosy) + (sin B — sin B' cos* y) (18) 
3 R / eo alia 
Fay = ida (y — y’ cos B) + 7 (sin y — sin y’ cos? f). (19) 


Equations (14), (15) and (16) can be surmised from Yamauti’s equations 
for vertical windows. One can also understand how equation (19) can be 
readily obtained from equation (18) by merely interchanging 8 with y and 
B’ with y’. 


In the equations derived above the angles are obviously to be expressed 
in radian measure. It however facilitates calculations to express the angles 
in degrees. Moreover it is usual to express the sky components as percen- 
tages. The equations modified accordingly are given below for sky com- 
ponents with the C.I.E. standard overcast sky. 


Vertical windows: 


_ B—B' cosy , sin (sin B sin y) 
*¥va = — 8-4 + 6-3 





— 4-547 sin 2y sin B’ (20) 


Fv. = ree + 9-095 (1 — cos B — cos y + cos 8 cos y’) 
(21) 


Fy, =P SY TY SINE + 9.095 (sin B— sin f’ cost y). (22) 





Horizontal windows: 





a f vty sin B 14.286 + 4-547 (sin 26 sin y’ 


+ sin 2y sin B’) 


__ sinc (sin 2 cos a_i (cos a cos y) (23) 
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Fug = =F + 9-095 (sin 8 — sin P’ cos? y) (24) 
'%F ay =F — B 9.095 (sin y — sin y’ cos? f). (25) 


[The angles in equations (20) to (25) are expressed in degrees.] 


The first terms in the above equations, with 3-6 replacing 8-4, gives 
the corresponding sky components with a uniform sky. 


The above equations for a point P at infinity get reduced each to a value 
zero. When the point P coincides with A, i.e., when 


B=y=B+y¥=5 or 90° 


the horizontal components get reduced each to a value 25%. The compo- 
nents in the vertical planes get reduced each to a value 


bh @ecy n 
(55+ 7) 100% = 19-8% (about). 


This later value can be easily derived as the sky component over either 
of the two mutually perpendicular vertical planes for a quarter hemisphere 
of a standard sky enclosed between them. For a uniform sky all values 
will be 25% when P coincides with A. 


The values calculated from these equations tally approximately with those 
given in Rivero’s tables, indicating thereby the care and accuracy with 
which the latter values have been computed though employing an indirect 
method. 


The equations given above are derived for a point P on the perpendi- 
cular through one corner of the window. If the point is not so situated 
and when the sill is above the point, the rectangular window should first 
be considered extended up to the foot of the perpendicular from the point 
on the plane of the window. The additional area can be considered as an 
algebraic sum of three rectangular areas each with one of their corners at 
the foot of the perpendicular from the point. The total sky component 
for the additional area, which can now be calculated, is then deducted from 
the component for the extended window, to give the component for the 
actual window. When the sill of the window is below the horizontal plane 
through the point, the portion of the window below that plane is neglected. 
It is thus possible to enlarge the utility of the equations to the general case 
of any point inside. 
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GLAZING CORRECTION 


The loss in transmission through glass at various angles of incidence 
can be calculated from the Fresnel’s equation. The loss is about 4% at 
normal incidence from each surface of the pane. At other angles of incidence 
the loss is more. Rivero (1958) has given the following empirical formula 
for the transmission factors at different angles of incidence. 


t™=c(l+sin* @)cos 6, where c=0-81 (about). 


This formula agrees with the theoretical curve fairly well, and is of a form 
convenient to adopt as an additional factor in the equations to be integrated. 
Sky components for glazed windows can then be easily derived using the 
general method adopted in this paper. Such equations are however not 
given here, as they are not likely to be of much practical significance. The 
glazing losses increase rapidly with the degree of dirtiness of the glass panes, 
depending on the period they are left uncleaned, the density of atmo- 
spheric dust at the locality, the inelination and orientation of the windows, etc. 
Such losses can be estimated only approximately and may amount to as 
much as 50% and more under bad conditions. There does not appear 
therefore to be any justification to calculate it with precision taking into 
consideration only the relatively minor variation of transmission factor 
with angle of incidence. Quite satisfactory results would be obtained for 
all practical purposes if the sky components for unglazed openings are 
reduced by a value 10 to 20% to allow for glazing losses, provided the glass 
panes are tolerably clean. 


SUMMARY 


Rigorous trigonometrical equations for calculating directly the sky 
components at any point due to a C.I.E. Standard Overcast Sky for both 
vertical and horizontal rectangular unglazed windows are derived and pre- 
sented. The equations, when suitably combined, can also be used for any 
sloping window. The general method adopted can be used to derive similar 
equations for a sky with any intensity distribution, provided that intensity 
is a function of altitude alone and independent of azimuth or vice-versa. 
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ABSTRACT 


The energy spectra of protons emitted in neutron bombardment 
of Al?’ is analysed in terms of the statistical model and volume direct 
interaction. It is found that the diffuseness of the nuclear potential or 
the form of level density used does not alter the shape of the energy 
spectra very much. The available experimental data agree fairly well 
with theory. The contribution of direct interaction is seen to be small 
and hence does not alter the general shape of the spectra. The direct 
interaction cross-section obtained from the measurements of the angular 
distribution in the case of iron and copper agree reasonably well with 
theory. 


INTRODUCTION 


In recent years a large number of experiments! *» * have indicated the 
inadequacy of the statistical theory in explaining the data. The two assump- 
tions of the statistical theory are (i) an intermediary compound system is 
formed which decays after a relatively long time, (ii) the contribution from 
the several overlapping levels in the highly excited compound system have 
random phases so that interference effects can be neglected. These assump- 
tions lead to symmetry about 90° in the angular distribution and linear 
dependence of the spectra of the emitted particles to the level density of the 
compound or the residual nucleus. The angular distribution would further 
be isotropic if the density of levels of total angular momentum J follows the 
(23 + 1) law.* 


Most of the experimental results indicate that (i) the angular distribu- 
tion is not symmetric about 90° and is mainly peaked in the forward angles. 
This anisotropy increases with increasing energy of the emitted particle, 
corresponding to the low-lying states of the residual nucleus. (ii) The high 
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energy region of the spectra is much larger than predicted by the level density, 
Also the spectra exhibit peaks in the high energy region which correspond 
to the discrete low-lying levels of the residual nucleus. There are some 
indications of the preponderance of low energy protons in (n, p) reactions.5 


The angular distribution of the emitted particles which leave the 
residual nucleus in discrete low-lying levels has been explained by the surface 
direct interaction which closely resembles stripping and pick-up reactions.* 
Unfortunately the theory is at present unable to predict absolute cross- 
sections. 


The general increase in the high energy part of the spectra may be due 
to any one of the following reasons: (i) Modification of level density 
formula taking into account of the effective mass of nucleons inside the 
nucleus as shown by the optical model.” (ii) Volume direct interaction.*"°. 
(iii) Interactions in which a large number of nucleons take part in the reac- 
tion like spot heating of the nucleus and collective excitation. 


The preponderance of low energy protons in (n, p) reactions may be 
due to the diffuseness of the nuclear potential," modification of the level 
density formula taking into account the effect due to pairing energy as sug- 
gested by Newton? and contributions from reactions like (m, np) due to 
compound as well as direct processes. 


In this paper we shall investigate the effects of these various factors on 
the shape of the energy spectra of protons emitted from Al®’ nucleus bom- 
barded with neutrons. It has been pointed out!* that volume direct inter- 
action does not contribute appreciably to the total direct interaction cross- 
section, especially in heavy nuclei, due to surface reflection and attenuation 
in nuclear volume. The effects of surface reflection and attenuation would 
be predominant if RA where R is the nuclear radius and A is the mean 
free path of nucleons in nuclear matter. Since A ~ 3-3 x 10-% cm. for the 
energy region under consideration,®»"* the above effects would not appre- 
ciably diminish the volume direct interaction cross-section for nuclei with 
mass number A ~ 60. To test this assumption, the experimental value 
of the direct interaction cross-section for iron and copper has been compared 
with the theory of volume direct interaction. 


THEORY 
Neutron bombardment of Al?’ would yield the following reactions: 
1. (n, p) and (n, pn) due to compound process. 
2. (n, n’), (n, n’ n) and (n, n’ p) due to compound process, 
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3. (n, p) and (n, pn) due to direct process. In the case of (n, pn) the 
neutron can be directly ejected simultaneously with the proton or can boil 
off the residual nucleus. 


4. (n, n’), (n, n’ n) and (n, n’ p) due to direct processes. In the case 
of (n, n’ n) and (n, n’ p) the second neutron or proton can be simultaneously 
ejected with the first neutron or can boil off the residual nucleus. 


5. (n, a) and (n, an) reactions. 


Reactions like (n, y), (n, d), (n, t), (n, 2p) are neglected either because 
their cross-sections for energies of incident. neutrons up to 20 MeV are 
negligible'® or because they are not energetically possible. 


The spectra of protons from (n, p) and (n, pn) reactions are given in 
the statistical theory by Blatt and Weisskopf.!® This expression is valid 
only when the level density follows the (2J+1) law. 


01 €0¢ (€) w (eg — €) de 


>F, (1) 





n(e)de = 


where o, is the cross-section for the formation of the compound nucleus 
due to the interaction of the target nucleus with incident neutron, «a¢ (¢) 
is the cross-section for the formation of the compound nucleus due to the 
interaction of the residual nucleus in the excited state with the emitted pro- 
ton of energy «, €, is the maximum possible energy for the emitted proton 
and involves the incident neutron energy and Q value of the reaction. 
w(e*) is the level density of the residual nucleus at the excitation energy, 
e*=e€,— ¢, D'F, is the total probability for the compound system 


to decay into the various channels v. Since (n, pp) for the energies under 
consideration is very improbable due to Coulomb barrier, we have assumed 
that only a neutron can boil off the residual nucleus. For the energies of 
interest we have assumed that simultaneous evaporation of two or more 
particles from the compound nucleus is improbable. As will be seen later 
the experimental data seem to support this approximation. 


The various cross-sections have been evaluated using the following two 
formule for the level density: 


w (e) = C exp. [2ate*] (2) 
w (€) = Cexp. [2a?(« — 4e)*] (3) 


where ‘a’ is a measure of the nuclear temperature and we have assumed it 
to be the same for all nuclei ander consideration. ‘C’ is a constant which 
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may vary widely from nucleus to nucleus and in effect it determines the 
level density for the particular nucleus. Je is the pairing energy for the 
nucleus under consideration and has been evaluated by Cameron.” 


Since it is very difficult to estimate the value of the constant C due to 
our lack of knowledge of level density in nuclei, we have eliminated it from 
the equations by using the experimental value of the total (n, p) cross-section. 


a;C . £ €0¢ (€) exp. (2ate’#) de 


[o (n,  — mar on = F, (4) 





where By is the binding energy of the last neutron in the residual nucleus 
and «’ = («,— «) or «’ = (e€,— « — 4e) depending on the level density 
formula used. [o(”, P)] comp. is the total cross-section for the (n, p) reaction 
going through the compound process. From (1) and (4) we obtain 

n (e) de = [o (, Ph, <%e (¢) exp. [22*<’#] de (5) 


J «oe (e) exp [2a*e’#] de 
€,—Ba 





whieh is independant of the constant C. [o(”, p)]comp. is obtained from 
the relation 


[6 , P) comp. = [7 PY ex. — [2 Hs PM a. (6) 


where [o (n, p)],,,, is the experimental value for the total (n, p) cross-section and 
[o (”, P)Jar. is the total (nm, p) cross-section due to direct interaction and is 
evaluated later [see eq. (14)]. 


The spectra of protons due to (n, np) reaction proceeding through the 
compound system is obtained as follows: 


Let N (en) den be the spectra of neutrons from the (n,n’) reaction. Let 
Bn and By be the binding energies of the last neutron and proton respectively 
of the residual nucleus. No further proton emission is possible if €, > 
(eon — Bp) where ¢gn is the maximum energy with which neutrons can be 
emitted in the (n, n’) reaction. Similarly no neutron emission from the 
residual nucleus is possible if en > (€9n — Bn). Let N (en, en’) denden’ be 
the number of secondary neutrons emitted in the range en’ and en’ + dey’ 
after the emission of N(en,)de, primary neutrons. Similarly let 
N (en, €p) dendep be the distribution function for secondary protons. From 
the statistical theory we have 


N (en, €n’) denden’ = Cn’ K (€n) o¢ (en’) exp. [2ateq”*) den’ (7) 
N (€n, €p) dendep = Cp K (€n) o¢ (€p) exp. (2a*en”*) dep (8) 
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where Cy,’ and Cy are the C-values of the corresponding residual nucleus 
after secondary emission. K(e,) is a constant of proportionality. ¢€,” 
and €p” denote the excitation energies of the corresponding residual nuclei. 
Since we assume that no other mode of decay is possible [the (n, n’ a) is 
neglected since the alpha particle, which has a relatively low energy, would 
have to overcome a much higher Coulomb barrier], 


N (€n) den = Cn'K (en) 7 at (€n’) exp. [2a*en”#] den’ 
0 


+ CpK (en) f epoe (ep) exp. (2aten" dep. (90 
° 





Hence 
K (€,)= 
N (€n) den 
Gt aim (en’) exp. [2a*en”*] den’ + Cp wy Zz Epc (€p) exp. [2atep”*#] dep 
(10) 
where 
En maz. = €on — €n — Bra. 
€p max. — fon — €n — Bp. (11) 


Hence if N (€n) den and Cy’/Cp are known, the (n, np) spectrum is obtained 
from equation (8). The total (m, 2m) cross-section in Al?’ cannot easily be 
determined since the ground state of Al®* has a very long half-life, (8-105y).% 
Hence we cannot eliminate Cy’/Cp from the above equation. For Al*’, 
the (n, 2n) reaction yields Al®* which is an odd-odd nucleus while the (n, np) 
reaction yields the even-even nucleus Mg**. Following the prescription 
suggested by Brown and Muirhead,’® we have assumed Cy’/Cp = 10. This 
assumption will not alter the spectra very much since the (n, 2n) contribu- 
tion in the energy region of interest is small, the threshold of the reaction 
being around 13 MeV. 


The (n, n’) spectrum can be obtained in the same way as the (n, p) 
spectrum and we get 





N (€n) den = [> a) EnFc (€n) exp. (2a*en”#) den (12) 


€n%c (€n) exp. [2a*en”*] den 


Where the lower limit of the integral is taken as zero since we include all 
reactions of the type (n, n’x). 
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[o (n, WT etic se [o (n, Pp) + o (n, PN) comp. 
— [o(n, p) + ¢ (a, pn)lax. — [o (1, 2) + ¢ (a, an)],,,, 
— [¢@,2')arJ (13) 


where one is the experimentally measured non-elastic cross-section. [o (n, p) 
+N, PN)]eomp, iS obtained from equation (5). [o(n,p) + (n, pn)],,,. and 
[o (2, n’)];,. from equations (14)-(16). [o(#, a) + o(”, an)],,.,,, has been 
evaluated from the experimental value of o(n,a) following the same pro- 
cedure as in the case of the (n,p)+(n, pn) reaction. We have chosen 
the total non-elastic cross-section since the (n,n’) cross-section is not known 
for these energies. Thus from (10), (11), (12) and (8) we can evaluate the 
proton spectrum due to (n, np) reaction. 


The contribution of the direct interaction process to these various 
reactions has been evaluated as follows. The nucleus is assumed to be 
composed of two non-interacting Fermi gases and the collisions of the inci- 
dent neutrons with individual nucleons in the nucleus give rise to these 
reactions. These collisions are considered to be free nucleon-nucleon 
collisions except for the restriction imposed on their final momentum states 
by the Pauli exclusion principle. We shall not give below the complete 
derivation of the various cross-sections but only quote the final results since 
these have been derived in detail by Hayakawa etal.® and by Brown and 
Muirhead.!° The neutron-nucleon collision inside the nucleus may yield 
any one of the following reactions: 


(1) »—~p collisions with both particles simultaneously ejected. 


(2) n —p collisions in which only the protons are emitted. 


In case (1) the final nucleus may boil off a neutron, thus yielding the 
(n,2np) reaction. For the energy region under consideration this is not 
energetically possible. In case (2) the final nucleus may similarly boil off 
a neutron yielding the (n, pn) reaction. Since we disregard the (n, pp) reac- 
tion, calculation of the (n, p) would automatically yield the (n, pn) reaction 
also. 


(3) n — p collisions with only the neutron emitted. This would yield 
a part of the (n,n’) process. The final nucleus may boil off a proton to 
yield the (n, n'p) or a neutron to give the (n, n’n) reaction. 


(4) n —n collisions in which either or both of the neutrons are emitted. 
The latter would give rise to (m,2n) reaction while in the former case the 
residual nucleus can boil off either a proton or a neutron to yield the 
(n,n’ p) or the (n, 2n) reaction respectively. 








AV id 
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If we obtain the effective (n,n’) spectrum due to (3) and (4) with one 
of the particles remaining inside the nucleus, then we can calculate the 
spectrum of protons due to (n,n'p) reaction with the proton boiling off the 
residual nucleus as in the case of the compound process [cf. eq. (10)]. 


These various cross-sections are given by: 











do 
iz) T (€p) ¢ (€p) 
—_ 3 
e- ex? [ Pp2np + Pn@nn x m 
i) _* (i) T (en) ¢ (en) [1 — T (ep) $ (ep)] (15) 
den/np i €,? [ Pomp + Pntnn x" 
+) -- (ane 2) Teas) — TP 


Xnp 
Fr" ce” | Pp P Xan P anp + Pnann | 


T(E) is the penetrability and ¢(¢) takes into account loss due to attenua- 
tion in nuclear matter and due to reflection at the nuclear surface. 
e-0-75R/Ae,) 


ae $(<) = [l—fi-T (2)} e-1- BRN) (17) 





where A(¢,) is the mean free path of nucleons inside the nucleus with energy 
(E+ V), V being the depth of the nuclear potential and E the energy of 
the nucleon outside the nuclear field. 


Xnp (E:) and Xnn(E,) are the total cross-sections for free n — p and 
n—n collisions with incident neutron energy E,. ’ 


All other notations are the same as in ref. 10. 


The functions (do/dx,'), (do/dx,"), (do/dx,'"") which occur in equations 
(14)-(16) are given by: 
- di . 
~ ) rim f F (Xp, X1) dx2 (18) 
xq” 
where x,” is x2’, x2’, OF xq/" are given by 
Xe = 1—(%,+Q — xp) 


X_° =. 1 — (4, — Xn) 
te = 1 — (x,’ —_ Xn) (19) 


~ 


A2 
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Q. 


Q —* 


x —= En + V ° xX y = En + V . 
— €p? 5] “nN ex” ’ 


er? = V — By; er” = V — Bn. 
By and By are the binding energies of the last neutron and proton in the 
nucleus respectively. 
3 1 -1 (2xot /2 (x1 + Xe) 
F (Xe, = = : tanh-! (= \A1 2 
ee ae x, + 3%, (21) 
In (21), x, is to be replaced by x,’ for the case x,=x,"’. Figure 1 gives 
values of (do/dx,”) for various values of x, as a function of x. 





We can define a general function a(x,) given by 


«(x)= f (aap) ave 2) 


where, x3 = 2 — (x, + Q’) for (n, p) and (2x — x,’) for (n, n’'). From equation 
(22) ann and ay, can be evaluated. Figure 2 gives a (x3) as a function of x, 
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for various values of x,. From these two figures the various direct inter- 
action cross-sections can be evaluated for any nucleus. 


Neutron and Proton Penetrabilities——The cross-section o,(<«) for the 
formation of the compound nucleus due to an incident particle of energy 
E with a target nucleus is given by 


o¢ (E) = 7k? z (2/ + 1) T, (EB) (23) 


where T; (E) are the partial penetrabilities as defined by Blatt and Weisskopf.'® 
The penetrabilities T(E) which occur in the above derivation for the proton 
spectra are defined as 


ee (E) 
TE) eR +A ate 
where A = (A/27), A being the wavelength of the incident particle and R is 


the incident channel radius. We have chosen for R the value of 1°5 A 
Fermis. 


For square well potential for protons o¢ (E) has been calculated from 
the computations of T;(E) by Feshbach et ail.'* 


The penetrabilities for charged particles have been calculated by Kikuchi! 
using the W.K.B. approximation. We have chosen the same potential 
for evaluating T(E). The potential is given by 

V (r) = Ve (r) + Vn (7) 
ZZ’e? 
r 


V,(r)=—V (* + xP. eS 2)" (25) 


where V, (r) and V(r) are the Coulomb and nuclear potentials respectively. 
We have chosen for the constants the values: a=41f; b=0°Sf. 
Then T(E) is given in the W.K.B. approximation by 


T(E) = exp. |- 2 Se f J/V®—E ‘| (26) 


Ve (r) = 





where r, and r, are obtained from the zeros of the function [V (r) — E]. 
Since the above approximation does not take into account the centrifugal 
force, it would yield a much higher value for the penetrability of high energy 
protons than given by the exact solution. To some extent this has been 
corrected for by smooth extrapolation of the above curve for low energy 
protons with the asymptotic value of T(E) given by 
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[T a ={1- aw. rl (27) 


(R +a) 


where B is the Coulomb barrier height for the nucleus. Figure 3 gives the 
penetrability as a function of proton energy for Mg?’ for square and diffuse 
potentials. 


RESULTS AND DISCUSSION 


The above theory has been applied to compute the energy spectra of 
protons emitted from Al®’ nucleus due to neutron interaction. The incident 
neutron energies chosen are 14-7 MeV and 17-5 MeV. The constants used 
in the computation are given in Table I. The value of ‘a’ which occurs 
in the level density formula has been varied from 1-5 Mev-! to 3 MeV 
in steps of 0-5 MeV-. 


Figure 4 gives a typical spectrum computed for an incident energy of 
14-7 MeV together with contributions from the various types of reactions 
considered above. It is seen from the figure that the contribution due to 
direct interaction is negligibly small and does not alter the shape of the total 
spectra appreciably except in the very high energy part. Since the differ- 
ential cross-section for high energy proton emission is small, experimental 
data have large statistical errors and hence cannot resolve the effect due to 
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direct interactions. The only other method of evaluating the direct inter- 
action cross-section from experimental data is to study the angular distri- 
bution. Allen? uses this procedure and assumes that the total non-iso- 
tropic part of the distribution is due to direct interaction. It has been shown 
by Ericson and Strutinski* that the angular distribution in the statistical 
model need not be isotropic and a more realistic assumption would be that 
it is of the form (A +Bcos?6). The experimental angular distributions 
are not accurate enough to permit such an analysis. 


Figures 5-8 give loge w(E) as a function of E, for the two types of 
potentials. and level density formulae that we have used. Here w (E) 
and E are defined as: 

N (€p) dep 


€p [ee ( €p)] square well 





w (E) = (28) 
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and E = Enp® — €p where E is the residual excitation energy if all the pro- 
tons were emitted through the (”,p) and (n, pn) reactions only. Enpp°® is 
the maximum proton energy possible for these reactions and ¢p is the energy 
of the emitted proton. N (€p) dep is the spectrum of protons emitted due 
to all the various types of reactions considered above. The reason for this 
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TABLE I 


Binding energy” of protons in: 











PP vctacvws isan neben cocce «65 BV 
5 Ae 8-3 MeV 
of neutrons in: 

Pe tinnbaaciaananiaes 7*7 MeV. 

De dikektiendtavien ean 13-0 MeV. 

DT ieakaieneth saben 6:4 MeV. 

Cross-sections 14-7 MeV. 17-5 MeV. 
a (n, p)5 53 mb.* 38 mb.* 
a (n, a) 120 mb. 80 mb. 
G5 960 mb. 900 mb. 
0, 860 mb. 830 mb.f 
Paring Energy”” 

SE” lub tirtiasinsetgiahe usbinarsb oie 2-1 MeV 

BT Achesitxcsecbadwendennia 4-3 MeV 

Se “sdkessaacbnebeernewvarines 2-2 MeV 


V = Nuclear potential = 42 MeV{ 


* This value corresponds to 60 mb. at 14 MeV and is lower than the one used in reference 1. 
This lower value seems to agree better with experiment and is the value obtained by Paul and 
Clarke™ and by Allen.’ 


t+ The o, has been calculated from the formula ¢, = #:(R + x)*. 

¢ This is the value obtained from the real part of the optical potential and seems to fit the 
total as well as the elastic neutron cross-section. Proton scattering data also yield a similar 
value. 
choice of definition is because the experimental data of Colli et al.5 is also 
represented in the same way. It is seen from Figs. 5-8 that the experimental 
points agree better for the case a=3MeV-! than for a= 1-5 MeV". 
The function loge w(E) is very insensitive to the type of potential used or 
to the level density formula assumed. The experimental points seem 10 
be in better agreement with the square well potential for 4E = 2-1 MeV 
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but no definite conclusion can be reached since the experimental errors are 
fairly large. The discrepancy between the 14:7 MeV data as pointed out 
by Colli etal.5 is only due to the (n, mp) reaction. 


The total cross-section for the proton production is more sensitive to 
the potential used and to the level density formula. Unfortunately we are 
unable to compare our results with the experimental data because of the 
large error introduced in the computation of the (n, 7’) cross-section from 
the total non-elastic cross-section. This can be avoided by a direct measure- 
ment of the (n,n’) cross-section. There seems to be some discrepancy in 
the total cross-sections for the various reactions but unless more accurate 
cross-section measurements are available nothing more conclusive can be 
derived from the present data. The deviation of the experimental data 
from the theoretical curves in Figs. 5-8 for large value of E may be due to 
errors in our estimation of the (n,n’) cross-section and our assumption of 
constant ‘a’ for all these nuclei. 


The experimental data for the differential cross-section for direct inter- 
action obtained by Allen! from angular distribution is compared with the 
theory in Figs. 9 and 10. The experimental points have been evaluated 
on the assumption that the statistical model yields isotropic angular 
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distribution. Since the data have large experimental errors, it is not possible 
to analyse it in terms of the A + Bcos?@ distribution. In the theoretical 
curve we have used the value of 0-25 for Xnn/Xnp, the ratio of free n —n 
to n — p cross-section. The actual value for this ratio lies between 0-25 
and 0-50!° and a higher value for this ratio would give a better fit to the 
experimental data. It seems from these figures that at least for nuclei in 
this region the major contribution to direct interaction is from volume 
effects. 


The above analysis indicates that a combined measurement of the 
various cross-sections as well as the energy and angular distribution of the 
emitted protons is needed for detailed analysis of these reactions. Also 
the present data are not accurate enough to decide between the types of 
potential or level density formula. 
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INTRODUCTION 


THE flow of heat through a roof slab is periodic, on account of the solar 
insolation, which imposes a diurnal variation of temperature on its exposed 
surface. The exact periodicity can be vitiated by the small difference that 
may occur between the initial and final temperatures for the 24-hour interval: 
but this difference is usually negligible, especially when the day forms part 
of a stable synoptic regime. The significant thermal property of the roof 
material: is not so much its thermal conductivity (K), as its thermal diffusi- 
vity, a given by the relation 


where p and Cp are its density and specific heat respectively. 


In common practice, a is calculated from the experimental values of K, 
p, and Cp measured under steady state conditions for oven-dry specimens. 
The temperature dependence of these constants may be a second order effect 
only ; but it is definitely known that they vary with the moisture content of the 
material. Since the air in the pores of a moist specimen is replaced by the 
more conducting water both in its liquid and vapour phases, not only the 
effective p, and Cp but K also alters. A vapour pressure gradient is set up 
due to temperature differences; and the consequent vapour movement 
inside gives rise to an additional heat transport, both latent and sensible, 
in parallel to the conduction flow. This reduces the thermal resistance or 
increases the apparent conductivity of the material. 


Considerable errors are therefore likely to be introduced in the analysis 
of heat transfer through building sections if a is calculated from the con- 
ventionally determined values of K, p and Cy. The correct approach would 
be to evolve a method of determining a of the material directly under the 
existing conditions of moisture. 
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An attempt is made in this paper to obtain a value for a for 7. roof slab 
from in situ temperature measurements of its outside and inside surfaces and 
of the ambient air inside the room below. 


EXPERIMENTAL 


Temperature data were available already, for the living rooms of about 
40 dwelling houses put up for the Low Cost Housing Exhibition at Delhi 
(Raychaudhuri, 1957). Thermocouples were set up at a number of points 
in the experimental rooms and hourly temperatures were recorded continu- 
ously for a number of days, to an accuracy of + 0-1°F. 


The data used in this paper were that obtained for one of the rooms, 
kept completely closed throughout a clear day (22nd April 1955) in the 
middle of a stable stretch of weather conditions. The roof of the room 
was a 44” R.C.C. slab. The room was 10’ 6”x 12’ x 9’ 9” in dimensions and 
built with 9” brick walls in cement mortar, plastered on both sides. 


Curves for the observed surface temperatures of the roof (outside) the 
ceiling (inside) and the ambient air temperature at the centre of the room, 
on the chosen day are given in Fig. |. 
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Fic. 1. Time-temperature curves for roof ceiling and indoor ambient air, (a) As observed, 
(b) As synthesised from the Fourier components. 
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THEORETICAL 


Houghten (1932) and his co-workers have given a solution for the ceiling 
temperature, 9¢ consequent on a unidirectional transverse flow of heat from 
the roof surface where the temperature varies sinusodially, with an ampli- 
tude @,. 


The solution is given below in a slightly altered form 
Ic = Rye~*.!8, sin (wt — B,l + 4) (1) 


where / is the thickness of the roof slab; 


8, -/2. (o=7 = 15°/hr.) 


_2V2 


| nee Sa 
1 VIF BA +AH 


fi 
= tan’ 5; ; 
° 2+h 


, h 
fh a KB, > 
and h is the appropriate film transfer coefficient, equal to the ratio between 


the rate of heat transfer per unit area of the ceiling and the ceiling-to-air 
difference in temperature. 


h, K and the ambient air temperature for the room below are assumed 
constant throughout the period while arriving at the above solution. 


This solution can be extended to the case where the roof temperature 
has any periodic variation other than simple sinusodial. 9. will now be 
given by the expression 


6. = Je + E RnePal8, cos (nwt — on — Bal + yn) (2) 


n=1 


when the roof temperature is given by the Fourier expansion 


0, = Oy + 3 On cos (nwt — on) 


ne1 
where 
Bn = VnBy 3 
2/2 





R, = SSS = 3 
"V2 + 2fn t+ fri 
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f 
= —1 Jn 
Yn = tan Fi7 
and 
hi 
fn ~ a/n 


8, and O¢ are related to 6g, the assumed constant air temperature by 
the steady state equation 


ci. « -~— 
] (0 — A) = h(O¢ — 4q) 


(Be — 8a) 


PROCEDURE 


The observed experimental time temperature curves of the roof (6,), 
the ceiling (@¢) and the indoor ambient air (q) are first analysed into their 
Fourier components as follows: 


6, = 86:07 + 23-02 cos(wt — 214° 9’) 
6°85 cos(2wt — 27° 59’) 

+ 1-01 cos (3wt — 115° 50’) 

+ 0-50 cos (4wt — 130° 22’) 

eet nn salnnaete (4) 

O¢ = 84:33 + 14-95 cos (wt — 242° 35’) 

+ 3-66 cos(2wt— 70° 30’) 

+ 0:32 cos (3wt — 163° 50’) 

+ 0-25 cos (4wt — 139° 31’) 

Pe sssccaeapcouaas (5) 
~-O_ = 81-45 + 6°50 cos(wt — 259° 54’) 

+ 1-76 cos(2wt— 84° 48’) 

+ 0-71 cos(3wt— 56° 44’) 

+ 0:34 cos (4wt — 117° 17’) 

mae - dubs Sigavee (6) 


-{- 
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The curves, summing up the first five terms of the Fourier expansions 
for 4», @ and 9g, are also given in Fig. | along with the observed curves. 
The correspondence between the two justifies limiting the calculations to the 
fundamental and the first three harmonics only and neglecting as insignificant 
the succeeding terms of the Fourier series. 


[It is to be noticed that the ambient air-temperature is not constant as 
assumed in the solution (Equation 2). However, if it is taken as constant 
at its mean value, 0g during the entire period, the solution for the ceiling 
temperature (Equation 2) can be applied as a first approximation, and more 
so when 0q = 0q = 81:45° F., at the instant defined by wt = 175° 16’ (as 
obtained from Equation 6). 


The next step would be to solve for a. But due to the inherent difficulty 
in employing the solution directly, an indirect method is adopted. The values 
of 4¢, the ceiling temperatures are computed at the instant corresponding to 
wt = 175° 16’, for different values of a from Equation 2, and substituting for 
h/K the value 1-609 as calculated from Equation 3. The relationship between 
a and the computed values of 9 at the specified instant (wt = 175° 16’) is 
shown graphically in Fig. 2. The observed value of ¢ at the instant is 90-82 
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Fic. 2. Computed ceiling temperature (when wt.=175° 16’) for different values of thermal 
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(as obtained from Equation 5), which when interpolated in the curve of 
Fig. 2 gives the value for a as 0:0322 sq.ft. per hour. 
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The ceiling temperature for the full cycle computed for this value of a 
can now be expressed by the equation: 


Ie = 84-33 + 14-86 cos( wt — 241° 33’) 
+ 3-59 cos(2wt— 76° 50’) 
+ 0:44 cos (3wt — 180° 15’) 
+ 0-18 cos (4wt — 207° 34’) (7) 


The curve of Equation 7 is given in Fig. 3, along with the experimental 
curve for the ceiling as defined by Equation 5. The correspondence between 


the two curves is good and the two coincide as they should when 0 = 9q. 
The discrepancy, where it exists, should be ascribed to the fact that the air 
temperature was not constant throughout as was assumed. Also the value 
of h can vary with time, being greater during the later hours of the night 
and early morning on account of a convective regime set up in the room 
air, 9 being less than 4g then (vide Fig. 3, where the room air temperature 
curve is also given for easy reference). During the rest of the period a strati- 
fied regime is set up as 0, > Og at that time. The reflections on the ceiling 
temperatures due to variations in %q and A will be greater during nights but 
in opposing directions. This will displace the observed curve above the 
calculated curve at first with a possibility of the two crossing over and the 
former coming below the latter after some time. During the day and up to 
about midnight the observed curve will be throughout above the calculated 
curve, as the effects due to / and 6g are additive then: but the shift may now 
be smaller in magnitude than that during nights as the relative variations in 
h and 9q from their average values are less during this period. These infer- 
ences, necessarily qdalitative, are borne out by Fig. 3. 


To a first degree of approximation at least, the value of a for the roof 
can be taken as 0-032 sq.ft./hr. based on in situ temperature measurements 
and the method of calculation adopted. This value of a is the average value 
of the thermal diffusivity of the roof slab for the complete cycle of 24 hours 
under existing conditions of moisture. 


SUMMARY AND CONCLUSION 


An attempt is made in this paper to compute the value of a for the roof 
material from in situ temperature measurements of roof and ceiling surfaces 
and indoor ambient air below, using the Fourier Components of the periodic 
variations and the available classical methods of solution for the general 
differential equation of unidirectional heat flow through the structure. Even 
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though the boundary conditions assumed in the solution are not strictly identi- 
cal with the experimental ones the computed thermal diffusivity, under in 
situ conditions, can be considered accurate enough for practical purposes. 
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(Communicated by Dr. S. Bhagavantam, F.A.sc.) 


1. INTRODUCTION 


THE measurement of porosity of a rock involves the evaluation of its pore 
volume and bulk volume. Any method that will enable the reliable deter- 
mination of these quantities may be employed. The bulk volume can, be 
determined either by direct measurement of its geometric form or by dis- 
placement of an equivalent volume of suitable liquid by immersing the coated 
specimen. The measurement of pore volume is an important and difficult feature 
of this experiment and several methods have been employed for this purpose. 
Tickell and others (1933) devised a porosimeter employing the phenomenon 
of change of length of mercury column with volume change in a manometer. 
Fancher (1950) employed the familiar Russel volumeter to measure the 
porosity of consolidated samples by saturating the sample with a suitable 
fluid and then observing the displacement or loss of weight resulting from 
immersion. Hughes and Cooke (1953) measured the porosity of rocks by 
using a modified Boyles Law apparatus which involves the evacuation of 
air contained in a rock by a pump. Davis (1954) suggested an empirical 
formula for estimating the porosity of sedimentary rocks from the bulk 
density. He had also shown that such estimated values closely agree with 
(3% discrepancy) the experimental results. The authors in the present 
investigation have employed Tickell’s porosimeter with certain important 
modifications which make the method simpler and the results repeatable. 
The ease and accuracy with which the pore volume could be determined 
for a given rock and thereby the porosity calculated from its bulk volume 
is the essential feature of this method. 


2. EXPERIMENTAL TECHNIQUE 
Figure | shows the modified form of Tickell’s porosimeter. The chief 
modifications introduced are: 


(1) The height of the specimen chamber is increased to handle bigger 
samples. 
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(2) The constancy of the mercury level at the bottom is maintained 
by keeping a small dish in another bigger dish and thereby allowing for an 
overflow, instead of a shallow dish employed by Tickell, where the level 
of mercury is not likely to be constant. 
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The bulk volume of the specimen is determined by measuring iis 
geometric dimensions. Closing the lower stop-cock, mercury is poured 
in the vertical tube to the mark (A). By closing the upper stop-cock and 
releasing ihe lower stop-cock, mercury is allowed to fall into the dish and 
attain a steady state which is taken as initial reading (1,). Again mercury 
is refilled to the mark (A) and in fact this operation is to be done whenever 
a fresh reading is to be taken. The dry specimen whose bulk volume is 
already determined is introduced in the specimen chamber (S). By releasing 
lower stop-cock, mercury is now allowed to come to rest at a higher position 
(1,), since there is a reduction of air in (S) due to the grain volume of the 
sample. The difference between !, and I, will be the length of mercury 
column corresponding to the grain volume itself. Care has been taken to 
see that joints are airtight and the instrument ts free from air-bubbles. 
From a knowledge of I,-L,, the grain volume is then read off from an 
already prepared calibration graph (Fig. 2). From the knowledge of the 
bulk volume and grain volume, the pore volume is calculated from which 
the effective porosity (P%) is computed. 
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The instrument is calibrated by taking a number of uniform ball-bearing 
shots (0:723 c.c.). As these B.B. shots are assumed to be practically non- 
porous, the rise in mercury column is due to the displacement of equivalent 
volume of air as a result of insertion of specimens in the chamber (S). 
Figure 2 shows the graph connecting the volume of ball-bearing shots with 
the corresponding rise in mercury level (I,-I,) for the apparatus we have 
constructed. It can be seen from the graph that every | c.c. increase of 
volume will result in rise of mercury column to 0-5cm. approximately. 
Further, the variation is found to be strictly linear as should be expected. 


The rock specimens chosen here are of definite geometric shape usually 
of rectangular parallelepipeds. Such specimens will enable us to deter- 
mine the bulk volume easily without taking recourse tc wax coating and 
subsequent measurement by displacing a suitable liquid equivalent to its 
volume as done by other workers. The present method enables us to repeat 
our experiments and get the mean values while in other methods this be- 
comes an impossibility. Another significant feature is that the specimen 
does not come into contact with the liquid. Air, because of its nearly 
perfect expansibility and low viscosity, penetrates the small and poorly con- 
nected pore spaces and gives a more accurate measurement of the effective 
porosity. This method is not affected by the variation of atmospheric 
pressure as we are taking I, and I,, at the same time and temperature cor- 
rection is made whenever required. 


3. RESULTS AND DISCUSSION 


Employing the porosimeter described above and taking all possible 
precautions, measurements have been made in some limestones and sand- 
stones. Results are given in Table I. Here density (p) is in grams/c.c., 
surface area (S.A.) is in square cm., Bulk Volume (B) is in cubic cm., dif- 
ference in mercury level [,-I, is in cm., grain volume (G) is in cubic cm. 
as read off from the calibration graph. Porosity (P%) is given as a 
percentage. 

I, is 21-2 at room temperature 27° C. and atmospheric pressure 760 mm. 


The vaiues (Table 1) are mean values of readings obtained a number of 
limes. The values are repeatable and the accuracy is found to be fairly high. 
The empirical formula suggested by Davis and referred to earlier for 


estimating the porosity (¢) in sedimentary rocks from Bulk density (pp) in 
CGS. Units is 


} — 


pa _— 
1-654 


| x 100. 















TABLE I 


Porosity values 
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Specimen p S.A. B I,-I, G Py 
sq.cm. Cc.C. cm. c.c. 

Limestone— 

1 2:7 8-98 1-818 0:8 1-7 6-49 

2 2°7 13-144 3-027 1-4 2-95 5:15 

3 2°7 35°897 10°85 5-1 10°65 1°84 

4 2-7 38-918 11-48 5-4 11-3 1-57 
Sandstone— 

1 2-47 24-138 5-615 2-1 4-3 23°44 

2 2-47 30-09 7-413 2°8 5-8 22:4 

3 2-47 34-687 9-05 3-4 7-1 21°55 

4 2:47 45-425 12-72 4-9 10-25 19-42 










equal to about 3%. 


TABLE II 


Comparison of porosity values 


The expected or mean error in ¢ was shown by him to be small and 
The authors have calculated the porosity employing 
this formula and the results (a) are tabulated in Table II along with the 
porosity values measured by the conventional method (b). Mean values 
obtained in the present investigation (c) are given for comparison. 






















Specimen a b c 
Limestone oo 2°05 1-96 1-84 
Sandstone .. 20°19 19-20 19-42 










The results in Table II are in good agreement with each other. 


In order to study the effect of size and surface area, number of 
parallelepipeds of limestones and sandstones have been prepared (perfect 
rectangular shape) and the results obtained are given in Table I. 
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Figure 3 shows the change in porosity value with surface area of the 
specimen. The graph shows that the porosity approaches constancy with 
increase in surface area in relation to bulk volume. Thus it can be seen 
that best results can be obtained by taking a number of small samples to 
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Fic. 3. Variation of Porosity with Surface Area. 


fill up the chamber rather than by taking one big sample. This is what one 
should expect. Care should be taken to see that the sample surfaces are 


plane parallel in order to get good results. 
4. CONCLUSION 


A simple porosimeter has been designed with a view to handle large 
volumes of samples and determine the grain volume with the help of a cali- 
bration graph. Bulk volume is determined directly from the geometric 
shape of the samples. The results obtained are repeatable and are found 
to be in good agreement with those obtained by other methods. 
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PHOSPHORUS oxychloride has been found to be a versatile solvating liquid 
for many inorganic compounds. A large number of non-polar halides of 
the elements of III, [V and V groups of the periodic table react with phos- 
phorus oxychloride and form complexes which are stable and easily isolated. 
Similarly a large number of solvates of metal oxides with phosphorus oxy- 
chloride have been reported in literature. This solvating property of the 
oxychloride is attributed to its co-ordinating tendency through either 
oxygen! or chlorine?-* acting as donor atom. 


Earlier reports about the complex formed between phosphorus oxy: 
chloride and aluminium chloride indicated the complex to be of composition 
AICl,.POCI,,* but recent observations concerning its composition are con 
flicting. Gutman® claims that the complex is AlCI,.POCI, while Groeneveld 
and Zuur® report the complex to be AICl,;.2 POCI, prepared under similar 
conditions. The latter authors have postulated the possible existence of 
other complexes with the following compositions: AICl,;.5 POCI, and 
AICI, .6 POCIs. 


The preliminary investigations of the present authors indicated that 
the molar ratio of aluminium chloride and phosphorus oxychloride in the 


complex formed was 1:2. A detailed investigation was therefore under | 


taken to throw more light on the composition of the complex and also to 
study its physico-chemical properties, such as its cryoscopic behaviour, 
electrical conductivity and absorption spectra in solution. 


EXPERIMENTAL 
Phosphorus oxychloride 


Phosphorus oxychloride of C.P. quality (E. Merck) was purified by 
distilling in an all-glass apparatus in an atmosphere of pure and dry nitrogen. 
The fraction distilling between 105-106°C. (683 mm.) was collected and 


preserved in sealed ampules (about 2-5g. in each). The purity of the | 
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oxychloride was checked by the estimation of the phosphate and the chloride 
obtained after hydrolysing a known quantity of the sample with dilute 
alkali. The results indicated that phosphorus oxychloride was 99-9% pure. 


Aluminium chloride 


Anhydrous aluminium chloride was prepared by passing dry hydrogen 
chloride over heated aluminium filings (Merck A.R. quality) in a pyrex 
glass tube maintained at 250°C. The condensed aluminium chloride was 
resublimed into a series of bulbs in a current of dry hydrogen chloride. The 
bulbs each containing about 2-3 g. of the anhydrous chlorides were sealed 
off. A representative sample was analysed by hydrolysing the chloride in 
a dilute solution of alkali and estimating the aluminium as the hydroxy 
quinolate and the chlorine by Volhard’s method. The analytical results 
showed that the anhydrous chloride was 99-8% pure. 


|. Preparation of the complex 

The method followed for the preparation of the complex was essentially 
the same as described by earlier workers with a few minor improvements. 
The transfer of materials, etc., was carried out in a ‘dry’ box and the 
experiments were conducted in all-glass vessels. 

A known quantity of anhydrous aluminium chloride was taken in a 
dry flask fitted with a phosphorus pentoxide guard tube. A known excess 
of phosphorus oxychloride was introduced into the same vessel. Consi- 
derable heat was evolved during the mixing of the reagents. The anhydrous 
aluminium chloride dissolved slowly in the oxychloride at room temperature. 
The reaction could be completed within fifteen minutes by heating at 95° C. 
The flask was cooled to room temperature and the excess of oxychloride 
was then removed in two ways: (i) by passing dry air through the solution 
and (ii) by removing the vapours under vacuum. 

(i) Sample obtained by passing dry air.—A portion of the solution of 
aluminium chloride in phosphorus oxychloride as obtained above was trans- 
ferred to a dry glass bubbler assembled along with suitable traps and guard 
tubes. Dry air was slowly drawn through the solution (for 24-36 hours) 
till a free flowing powder was obtained. The free flowing powder samples 
(2-3 g.) were preserved in sealed tubes. 

(ii) Removal of the solvent by evacuation—Another portion of the 
homogeneous solution was transferred to an evaporating dish which was 
placed in a vacuum desiccator and evacuated slowly at 25°C. with occa- 
sional stirring till a free flowing powder was obtained. The powder samples 
were preserved in sealed tubes as mentioned above, 
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Analysis of the samples.—The estimation of aluminium and phosphorus 
in presence of each other presents some difficulties. In order to obtain 
reliable analytical results the method described in an earlier communication 
was adopted for the determination of both phosphorus and aluminium? 
The chloride content of the samples was determined after hydrolysing the 
substance with a cold dilute solution of sodium hydroxide by Volhard’s 
method. This procedure prevented any loss of hydrogen chloride during 
hydrolysis in water and gave correct analytical results.® 


The analytical results indicated that a stable complex is obtained in 
whatever way the excess of oxychloride is removed from a solution of alu- 
minium chloride in excess of phosphorus oxychloride. The molecular ratio 
of aluminium chloride and phosphorus oxychloride in the complex was 
found to correspond to 1:2. 


As has been pointed out earlier, Gutman reports a complex of the com- 
position AICI,.POCI,. But no details of the conditions of preparation and 
methods of removal of the excess oxychloride are described. Under the 
specified conditions described above, only one type of complex can be iso- 
lated corresponding to the composition AICl,.2 POCI;. The physico-chemical 
properties of the complex are described below. 


2. Cryoscopic behaviour of the complex in nitrobenzene 


Preliminary investigations indicated that nitrobenzene was the most 
suitable solvent for the cryoscopic work. The solubility of the complex 
in nitrobenzene at 25°C. was determined and this corresponded to the 
value 5-91 g. in 100g. of nitrobenzene® whereas other organic solvents like 
carbon tetrachloride, benzene, etc., could dissolve both aluminium chloride 
and phosphorus oxychloride separately but the complex was found to be 
insoluble in such solvents. Therefore all the cryoscopic measurements were 
carried out in nitrobenzene solutions only. 


Preparation of pure and dry nitrobenzene.—Nitrobenzene (C.P. quality) 
was purified by the method similar to that described by Meisenheimer and 
Dorner’® and Roberts -and Bury.“ A pure sample of nitrobenzene was 
distilled and the middle portion was fractionally crystallised three times, 
dried over anhydrous calcium chloride and again distilled. The middle 
portion was stored over phosphoric anhydride and just before commencing 
the experiment was distilled at 10mm. pressure into cryoscopic tubes. A 
carefully dried all-glass distillation apparatus was used for such purposes. 
The freezing point of such a sample was 5-77° C. + 0-005. 


The purity of nitrobenzene and its cryoscopic constant (K+) were checked 
by determining the molecular weight of resublimed naphthalene which is 
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known to be a monomer in this solvent." Ky of nitrobenzene was taken 
to be 6°89! which is in good agreement with the value calculated using 
latent heat of fusion of nitrobenzene as 22-46 cal./g. as determined by 


Meyer.” 


The molecular weights of anhydrous aluminium chloride, and phos- 
phorus oxychloride as well as that of the complex were determined separately 
and the following values were obtained: 


Aluminium chloride: Mol. weight in nitrobenzene 134-6 + | 
Phosphorus oxychloride = 163-7 +3°5 
“Complex” 345°6 + 5 


It has been observed that aluminium chloride has the formula -Al,Cl, 
up to 400°C. corresponding to all the three states of aggregation.'* The 
molecular weight of aluminium chloride in solutions of liquids like carbon 
disulphide, benzene, etc., also corresponds to the dimeric form; whereas 
in solvents like pyridine and ether! the anhydrous chloride is essentially 
monomeric. In general, it has been claimed that in any solvent that is 
known to combine with anhydrous aluminium chloride, aluminium chloride 
has a molecular weight corresponding to AICls. 


Aluminium chloride forms an addition product with nitrobenzene also'® 
with the formula (AICI,.Cs,H;NO,). The molecular weight of anhydrous 
aluminium chloride as determined in the present experiment indicates that 
it is in monomeric form in nitrobenzene solution at comparatively low 
concentrations. This result is in agreement with the values reported in 


literature.” 


Phosphorus oxychloride has a molecular weight in nitrobenzene corres- 
ponding to the monomeric form. While at higher concentrations there 
seems to be a tendency towards dimerisation. However, the values obtained 
in the present investigation agree better with the theoretical value than 
reported earlier.1® 


The determination of the molecular weight of the complex formed 
between phosphorus oxychloride and anhydrous aluminium chloride in 
nitrobenzene has been carried out using comparatively dilute solutions. In 
spite of this the value obtained in the present investigation (345-6) is far 
from the theoretical value corresponding to the molecular formula AICI. 
2 POCI, (440-1). This deviation of the experimental value from the theore- 
tical value may be attributed to the dissociation of the complex at low con- 
centration into phosphorus oxychloride and a lower complex in the following 
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way: AICI,.2 POCI,; = AICl,;.POCI,; + POCI;. One may even expect 
the complete breakdown of the complex into aluminium chloride and phos- 
phorus oxychloride thus giving rise to three molecules instead of two mole- 
cules as suggested in the foregoing. A comparison with the complexes 
formed between zirconium and hafnium tetrachlorides with phosphorus 
oxychloride eliminates the probability of getting three molecules.!® 
3. Electrical conductivity in nitrobenzene 

In order to ascertain whether the complex corresponding to the mole- 
cular composition of the solid existed in solution, electrical conductivity 
measurements of the complex in nitrobenzene solution were carried out. 
A modified form of Freas conductivity cell was used. The electrodes con- 
sisted of two circular platinum plates (diameter ?”) kept apart at a distance 
of 2mm. The conductivity cell had an all-glass airtight lid in order to 
protect the solution from contamination with atmospheric moisture. The 
measurements of conductivity were carried out with the help of ‘Electro. 
spezial” conductivity meter (Model GM 4249). The conductivity cell was 
kept in a thermostat at 25 + 0-05° C. 


The electrodes were platinised and the cell constant was determined 
in the usual manner. The cell constant had a value of 0-03865 which enabled 
the measurement of conductivity of the low conducting non-aqueous solutions. 
The nitrobenzene used had a specific conductivity of 0-6 x 10-7 mhos per 
cm. at 25°C. and this agrees with the value given by Payne (1 x 10-7? mho 
per cm.) % @ 


Aluminium chloride in nitrobenzene had a specific conductivity of 
4x 10> mhos/cm. at a concentration of 0-0! mole per litre which increased 
to 148-7x 10-5 mhos with increase in concentration to 1-01 mole per litre. 
Similarly the specific conductivity of phosphorus oxychloride solution in 
nitrobenzene increased from 1-69x10-® mhos/cm. to 25-43x10-* mhos 
when the concentration was increased from 0-02 mole/litre to 2-13 moles/ 
per litre. These results are presented graphically in the curves A and P of 
Fig. 1. These experimental values are in agreement with the results reported 
in literature.1§» 24 


The electrical conductivity of different mixtures of aluminium chloride 
and phosphorus oxychloride in nitrobenzene prepared by mixing the individual 
components to give suitable molar ratios at different concentrations were 
measured. The change in conductivity of the mixture of the reactants in 
nitrobenzene while keeping the concentration of aluminium chloride complex 
isolated (AICI,.2 POCI,) were also measured. 
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It is known that anhydrous aluminium chloride conducts electricity 
in all the three states of aggregation except at the melting-point. The ionic 
species responsible for conducting the current are reported?* to be Al,Cl, 
= Al*3 + (AICI,)-*. In solvents like ethyl chloride the aluminium cation 
is observed?* to be heavily solvated® as [Al(C.H;Cl), [**. It is also known 
that in aqueous medium the aluminium ion is solvated corresponding to a 
co-ordination number six. From the cryoscopic data, it has already been 
indicated that the monomeric form of aluminium chloride, is associated 
with nitrobenzene but the extent of solvation is not known. It is very 
difficult to distinguish from the cryoscopic data whether the two particles 
correspond to 2 AICI, or Al*® and (AICI,-*) owing to autoionisation. It 
is reasonable to attribute the electrical conductivity of aluminium chloride 
in nitrobenzene to these ionic species. 

Phosphorus oxychloride is not a very good conductor in nitrobenzene 
when compared with aluminium chloride. Even this feeble conductance 
which is not very much different from the solvent itself may be attributed 
to the ionic species formed by the autoionisation of phosphorus oxychloride 
into: 2 POCI,=POCI,* + POC], as suggested by earlier workers,’ 4 
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Many of the electrolytic properties of phosphorus oxychloride have been 
explained on the basis of such an ionisation. 


The conductivity of various mixtures of aluminium chloride and phos- 
phorus oxychloride in nitrobenzene in different proportions is higher than 
the sum of the individual conductivities of the two components calculated 
from the graph. For instance, the specific conductivity of a mixture of 
equal volumes of 0-6 M AICI, and 0:6 M POCI, (the resulting mixture will 
have 0-3M AICI, and 0-3 M POCI,) was found to be 11-04~x 10-* mhos/cm. 
which is nearly double the value of the sum of the individual conductivities 
at the corresponding concentrations in solution (6-75 x 10-* mhos for AICI, 
0-079 x 10-* for POCI, 6-83 x 10-4). 


In another case when equal volume of 0-4M solution of aluminium 
chloride was mixed with 0-6M solution of phosphorus oxychloride in 
nitrobenzene, the conductivity was found to be 12-47x10-* mhos/cm. 
Whereas the sum of the conductivities of individual components corres- 
ponding to 0-2M aluminium chloride and 0:3M phosphorus oxychloride 
is equal to 4-78 x 10-* mhos/cm. 


Similar observations are made in other cases also. It can be inferred 
from the above results that the contribution of phosphorus oxychloride to 
the electrical conductivity is very small when compared with that of alu- 
minium chloride. However, the mixture shows a marked increase in the 
conductivity indicating thereby the conducting species are very much 
different from either AIC], or POCIs. 


The conductivities of the solution (Table I) prepared from the isolated 
complex AICI,.2 POCI; in nitrobenzene also show a higher value than the 
sum of the conductivities of the individual components at corresponding 
concentrations. However, the solubility of the isolated complex,® which 
corresponds only to 0-53 molar, restricts the choice of higher concentration. 
Heating the mixtures of the two components in nitrobenzene solution to 
95° C. had little effect on conductivity indicating that the reaction between 
aluminium chloride and phosphorus oxychloride on mixing was complete 
even at room temperature. 


4. Conductometric titration of aluminium chloride with phosphorus oxychloride 
in nitrobenzene 


A solution of aluminium chloride in nitrobenzene (0-37 M) was titrated 
against a standard solution of phosphorus oxychloride in nitrobenzene con- 
ductometrically and it was found that the specific conductivity of the 
resulting system attained a limiting value at a point corresponding to the 
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TABLE I 
Electrical conductivity of the complex AICI,.2POCI, in nitrobenzene at 25°C. 





Concentration in mole per litre 0-34 0-68 2-05 3-42 4-79 6:16 6-84 
x 10? 


Sp. conductivity in mho./percm. 10* 0-01 0-27 0:98 1°53 2-04 2-69 3-18 








molecular ratio of the components in solution AICIl,: POCI, 1:2. Con- 
ductometric titration curve (C) is given in Fig. 1. This value clearly 
indicated the formation of the complex corresponding to a definite compo- 
sition AICl;.2 POCI,. 


It has already been pointed out that the conductivity of the complex 
in solution is very much higher than the sum of the conductivities of the 
two components at corresponding concentrations. This may be attributed 
to the formation of AICl,- and POCI,* ions and the tendency of the additional 
phosphorus oxychloride to dissociate into POCI,;* and Cl- ions. All these 
species are likely to bring about a perceptible increase in specific conductivity 
of the solutions of the complex formed between aluminium chloride and 
phosphorus oxychloride in nitrobenzene. 

5. Absorption spectra of the systems AICl;-POCI, in nitrobenzene 


The absorption spectrum of pure and dry nitrobenzene showed a peak 
at 4150 A with an optical density of 9-25x10-*. This is in agreement with 
the value reported in literature.2* Absolute alcohol was used as a standard 
while measuring the absorption of nitrobenzene. The absorption spectrum 
of phosphorus oxychloride in nitrobenzene, using nitrobenzene itself as a 
standard, did not show any characteristic absorption maximum. 


Pure and dry nitrobenzene has a straw yellow colour, whereas even a 
dilute solution of aluminium chloride is deep yellow. The intensity of the 
colour depends on the concentration. The absorption maximum is found 
to shift towards longer wavelengths with increase in concentration. The 
absorption maximum changed from 4400 A to 4620 A when the concentration 
of the solution was changed from 0-02 M to 0-4M. Beer’s law was found 
to be valid between concentrations 0-016M and 0-025 M of aluminium 
chloride in nitrobenzene at a wavelength of 4400A. The optical density 
of aluminium chloride solution in nitrobenzene reached a limiting value of 
0-825 at 0-22M. 


The optical density measurements of the mixtures containing various 
proportions of aluminium chloride and phosphorus oxychloride (corres- 











278 K. N. VENKATARAMAN AND A. R. VASUDEVA MURTHY 


ponding to solutions of concentration in the range where Beer’s law is valid) 
indicated that there was another range of maximum absorption between 
6950 A and 7200 A which is not observed in case of either nitrobenzene or 
the individual solutions of aluminium chloride or oxychloride. It could 
therefore be inferred that there was complex formation between aluminium 
chloride and phsophorus oxychloride when mixed in nitrobenzene solution 
even at room temperature. 


The solutions of the isolated solid complex in nitrobenzene also exhibited 
a similar peak in the range 6900 A-7200 A. The complex whether it is 
formed in solution by the addition of aluminium chloride and phosphorus 
oxychloride in nitrobenzene or prepared separately and then dissolved in 
nitrobenzene appear to be the same. 


SUMMARY 


1. Anhydrous aluminium chloride reacts with phosphorus oxychloride 
to give a complex with a composition AICl;.2 POCI, which can be prepared 
in the form of a free flowing powder. 


2. The phosphorus oxychloride-aluminium chloride complex in nitro- 
benzene dissociates into AlCl,;.POCI,; and POCI; as indicated by the cry- 
oscopic measurements. 


3. The solution of the complex in nitrobenzene has a higher specific 
conductivity than the corresponding electrical conductivities of individual 
components. Similar higher electrical conductance is observed when the 
two components are mixed in nitrobenzene in different proportions. 


4. When a solution of anhydrous aluminium chloride in nitrobenzene 
is titrated conductometrically against a solution of phosphorus oxychloride 
in nitrobenzene, a limiting value in the conductivity is reached at point 
corresponding to the molecular composition, the components in the ratio 
of 1:2 AICl,: POCI, in solution. 


5. The absorption maxima of the complex in nitrobenzene solution 
differ from the absorption maximum of the individual components. 
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INTRODUCTION 


THE study of clay minerals is of comparatively recent origin in the develop- 
ment of geological sciences. On account of the importance of those minerals 
in agriculture, construction engineering, and also in ceramic, rubber, paper 
and petroleum and other industries, considerable information has been 
gathered. Recent and modern research techniques like X-ray, Differential 
Thermal Analysis, Electron Microscopy and Infra-red Spectroscopy have 
opened new vistas for fundamental study. 


The clay materials available for study are mainly of two kinds—recent 
clays, unconsolidated or in the process of consolidation, and ancient clays 
from the sedimentary formations of the geological past. In the case of the 
former, the environment of deposition is readily known. The environment 
of ancient sedimentary clays can be ascertained by other petrographic methods. 
A critical study of the clay minerals with reference to their environment of 
formation, and the changes that are brought about by the later processes 
of diagenesis and lithification, will be interesting. The present investigation 
is an attempt at the study of clay minerals from their environmental point of 
view. In this paper, the results of Differential Thermal Analysis of eleven 
clay mineral extracts from shales and shaly sandstones belonging to the 
Chintalapudi formation, equivalent to the Kamthis of Triassic age, from 
the Lower Gondwanas of Krishna District and parts of West Godavari 
District, in Andhra Pradesh, are presented. 


LOCATION AND NATURE OF THE SAMPLES 


The Chintalapudi formation in the Krishna District occurs as a small 
patch of about five square miles near Somavaram (80° 55’ ; 16° 55’). Seven 
of the samples are from this patch. 
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The other four samples are from the outcrops of the formation east 
of Chintalapudi (81° 3’ ; 17° 3’). The outcrops in fact form part of the 
Kamthi beds of Nagpur and Chanda, extending down the Godavari valley. 
The Somavaram and the Chintalapudi areas are removed from each other 
by about 15 miles and separated by khondalites. The rocks, however, bear 
a close resemblance and are composed of colourless or light coloured, highly 
friable, medium to coarse-grained felspathic sandstones and fine well-lami- 
nated shales, with intercalations of hard reddish ferruginous streaks. A 
number of fossil impressions of Glossopteris and Vertebraria have been 
described from both the areas and determine the age of the beds. The older 
formations in the vicinity from which the Chintalapudi formation could 
have derived their material, are quartzo-felspathic, sillimanite nee and 
schists and quartzites of the khondalite suite. 


EXPERIMENTAL PROCEDURE 


All the samples were crushed and the coarse fraction separated by 
sieving. After dispersal in 0-01 N Sodium hexametaphosphate solution 
for a few days, a further separation was carried out by wet seiving using a 
phosphor-bronze sieve. Finer fractions were obtained by settling methods. 


The Differential Thermal Analysis apparatus used is of a portable type 
and the results are qualitative. The instrument consists of two junctions 
of a differential thermocouple made of chromel : alumel, placed symmetri- 
cally in two holes of a metal crucible block and connected to a sensitive 
galvanometer, and a pyrometer connection placed in the same block at right 
angles to the line joining the thermocouple junctions. The substance to 
be studied and calcined alumina, which is taken as the inert substance, are 
packed round the two junctions. The whole crucible is heated by an elec- 
trical furnace, its rate of rise in temperature being controlled by a variable 
resistance. Sufficient care is taken to see that uniformity in nature of packing 
and in rate of rise of temperature at 15°C. per minute are maintained. 


The instrument was initially calibrated with pure quartz powder. The 
difference between the ideal inversion temperature of quartz (573° C.) and 
the temperature at which inversion actually took place, under the existing 
conditions as recorded by the pyrometer of the instrument, was found to be 
20°C. and is applied as a correction to all the pyrometer readings. 


RESULTS 


As can be seen from the differential thermal curves shown in Figs. | and 
2, samples 1-7 show an initial endothermic peak between 75° C. and 150° C., 
another between 500° C, and 620°C, and an exothermic reaction at 920°C, 
A4 
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Samples Nos. 2, 4, 5 and 7 show, besides, double peaks in the second endo. 
thermic trough. Samples Nos. 8, 9, 10 and 11 are similar in their behaviour 
to the previous ones, except for the absence of the initial endothermic peak. 
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The curves obtained are compared with standard curves of different 
clay minerals given by Grim and Rowland (1942) and Kerr and Kulp (1948). 
The samples studied correspond very closely with minerals of the kaolinite 
group. The minor variations noted may be due to the occurrence of the 
clay minerals as mixtures and due to the presence of some impurities like 
sericite. 

DISCUSSION 


Frederickson (1952) points out that kaolinite characteristically develops 
in an acid environment. Méillot (1952), who has analytically studied some 
sedimentary clays with reference to their environment of deposition, comes 
to the conclusion “‘ Kaolinite sediments are due to the decomposition of 
crystalline rocks during the weathering. The rock must have been crystalline, 
because the waters that come from a mountain range covered by a sedimentary 
carapace are basic and are unfavourable to the formation of kaolinite, 
Carbonates are completely absent which is a useful criterion ”, 
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Grim, Dietz and Bradley (1949), based on their work on the Pacific 
sediments, point out that kaolinite is lost during marine diagenesis, but add 
that clay mineral investigations of recent and ancient sediments are too scanty 
to determine the general validity of this conclusion. However, since many 
definitely marine samples studied by these authors and others show illite 
and montmorillonite, rather than kaolinite, to be the dominant clay mineral 
component, an abundance of kaolinite may safely be taken to indicate a 
non-marine origin for the deposit containing it, in paleogeographic inter- 
pretations of an ancient sediment. 


Krumbein and Sloss (1953) state that kaolinitic clay minerals are- the 
normal end-products of podsolisation and lateritisation. There is also a 
tendency in this process of soil formation for leaching iron compounds. 


It has been mentioned earlier that the source rocks for the Chintalapudi 
formation appear to be the para-gneisses and schists of the khondalite suite 
of rocks. The author has observed that, in the areas under consideration, 
the felspar content of the gneisses and schists is comparatively more than in 
some other khondalite areas. Weathering of these would have yielded the 
kaolinitic group of clay minerals. The presence of plant fossils to the entire 
exclusion of marine fossils can also be taken as an evidence for continental 
conditions of deposition. The streaks of ferruginous matter occurring as 
intercalations in the Chintalapudi formation are probably the result of the 
natural tendency for leaching of the iron compounds under podsolisation. 


Size analysis of the Chintalapudi samples indicates variations in the 
sorting of sediments. There is better sorting towards the centre of the patch 
of the formation compared to the marginal areas. This is generally true 
in the case of lakes. 


CONCLUSIONS 


The Chintalapudi samples studied show the presence of the kaolinite 
group of clay minerals. Though the samples are from two different areas, 
there is a consistency in the clay mineral composition. It may be concluded 
that the formation was deposited under non-marine conditions, possibly 
in a lacustrine or fluvio-lacustrine environment. X-ray data, which will 
be presented elsewhere, confirms the identification of the minerals. 


SUMMARY 


Eleven clay mineral extracts from shales and shaly sandstones equivalent 
to Kamthis of Triassic age, from the Somavaram (80° 55’ ; 16° 55’) and 
Chintalapudi (81° 3’ ; 17° 3’) areas, are studied by Differential Thermal 
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Analysis. The presence of kaolinite group of clay minerals is indicated, | 
The kaolinites would have formed by the erosion of the para-gneisses and | 
schists, which are the source rocks for the Chintalapudi sandstones, and 
are suggestive of deposition under continental conditions, possibly in a 
fluvio-lacustrine environment. 
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